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Abstract 

We study the theory of systems with constraints from the point of view of 
the formal theory of partial differential equations. For finite- dimensional 
systems we show that the Dirac algorithm completes the equations of 
motion to an involutive system. We discuss the implications of this iden- 
tification for field theories and argue that the involution analysis is more 
general and flexible than the Dirac approach. We also derive intrinsic 
expressions for the number of degrees of freedom. 

1 Introduction 

Constrained dynamics represents a cornerstone of theoretical physics, as every 
relativistic theory and every theory with gauge symmetries necessarily possesses 
constraints. Thus it is not very surprising that many methods for dealing with 
such systems have been developed (see e.g. Jl^, |l5|, [l9[ ||[ Q). The purpose of 
this and the following articles in this series is to present an alternative ansatz 
based on the modern theory of differential equations and especially on the con- 
cept of involution |23|, ||^] . 

The classical Hamiltonian treatment of systems with constraints was devel- 
oped by Dirac ||, |j| |(J. We will show that in the case of finite-dimensional 
systems his algorithm corresponds to rendering the equations of motion involu- 
tive. In the language of exterior differential systems this was already noted by 
Hartley et al. \L^ . But we will also show that this connection does no longer 
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hold for field theories. Here it might happen that the Dirac analysis alone is 
not sufficient to obtain all constraints. 

This identification appears natural, as the basic idea behind the Dirac al- 
gorithm is to check whether or not the equations of motions are consistent. 
But the notion of involution represents essentially a mathematical formulation 
of this problem. We believe that this approach has conceptual and practical 
advantages, especially in the case of field theories. The theory of involution is 
well understood for arbitrary systems. The Cartan-Kuranishi Theorem [Q |3| 
yields a general procedure to complete any system of partial differential equa- 
tions to an involutive one. 

In contrast one can find in physics many different approaches depending on 
whether one deals with a system in Lagrangian or Hamiltonian formulation, 
or whether the Lagrangian contains higher order derivatives or is linear in the 
velocities. Each case is handled individually in the literature. A closer analysis 
shows, however, that almost all proposed methods represent nothing else than 
special cases of the general completion procedure stemming from the Cartan- 
Kuranishi Theorem^. 

We study in this article the standard situation of a system described by 
a Lagrangian depending on the generalized coordinates and the velocities as 
well as Lagrangians containing higher order derivatives. The special case of a 
Lagrangian being linear in the velocities and the so-called symplectic formalism 
of Faddeev and Jackiw || will be considered in the next article in this series. As 
concrete examples we consider among others the rigid rotator and Podolsky's 
generalized electrodynamics. 

Special emphasis is put on the problem of determining the number of degrees 
of freedom. We propose a new intrinsic definition for field theories using the 
Cartan characters of the field equations. It can also handle systems described 
in characteristic coordinates like light-cone coordinates. The classical approach 
based on a distinction into first and second class constraints fails in such a 
situation J32], as seemingly too many constraints occur. 

The article is organized as follows: The next two sections serve as a brief 
introduction into the formal theory of differential equation. They define the 
notion of involution and show how one completes an arbitrary system to an 
involutive one. Section || reviews the classical Dirac approach. In Section ^| we 
show its relation to the formal analysis of the Hamiltonian equations of motion 
for finite-dimensional systems. A detailed example is considered in Section ^. 
Section |^ contains an example of a field theory where the Dirac algorithm alone 
does not suffice to exhibit the full constraint structure. The problem of counting 
degrees of freedom is tackled in Sections || and |To| for the finite and infinite- 
dimensional case, respectively. Finally, before some conclusions are given, we 
consider in Section ^ Lagrangians depending on higher order derivatives. 

x This holds even for the approach of Lusanna [juj via the Noether theorem, as it is based 
on the eigenvectors of the Hessian. But these correspond exactly to the linear combinations 
of the equations of motion which yield the integrability conditions. 
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2 Involution 



Formal theory uses a geometric approach to differential equations based on the 
jet bundle formalism. It is beyond the scope of this paper to give a detailed 
introduction into the underlying theory. The interested reader is referred to the 
literature J23|, ^6|, ^8[ . Here we are concerned with two topics: the definition 
of an involutive system and how to compute the arbitrariness of the general 
solution of such a system. 

We will always work in a local coordinate system, although the whole theory 
can be expressed in a coordinate free way. Let X denote the space of the 
independent variables x±, . . . , x n and let the dependent variables u u m be 
fiber coordinates for the bundle £ over the base space X. Derivatives are written 
in multi-index notation p® — d'^'u 01 jdx^ ■ ■ ■ dx£ n where \fj,\ = (j>x + • ■ ■ + fj, n is 
the length of the multi-index \x = [hi, . . . , /j, n ]. Adding the derivatives p°^ up to 
order q defines a local coordinate system for the q-th order jet bundle J q £ . A 
system of differential equation !Z q of order q can be described locally by 



K 



{ 1> T (x. l ,u a ,p°)=0, r=l,...,p; \fx\<q. (1) 



Geometrically, this represents a fibered submanifold of J q £ . 

At least some of the ideas behind the concept of involution can be under- 
stood best by considering the order by order construction of a formal power 
series solution. For this purpose, we introduce the symbol M q of a differential 
equation lZ q . 

Definition 1 The symbol A4 q of the system ^) is the solution space of the 
following linear system of (algebraic!) equations in the unknowns 



M 9 : £ -5^K=0- ( 2 ) 




(By abuse of language, we will refer to both the linear system and its solution 
space as the symbol). 

The placeholders are coordinates of a finite-dimensional vector space, i.e. 
we introduce one coordinate for each derivative of order q. Definition |l| is most 
easily understood by considering a quasi-linear system, i.e. a system linear in 
the derivatives p^ with = q. For such a system the symbol is simply obtained 
by taking only the linear highest order part and substituting for p". 

We make a power series ansatz for the general solution of the differential 
equation lZ q by expanding around some point x° 

oo a 

E -$(*-*°r (3) 
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and substitute this ansatz into the equations (0) evaluating at x° . This yields 
a system of algebraic equations for the Taylor coefficients a" up to order q. 

The remaining coefficients can be computed by linear algebra only. For 
the coefficients of order q + r we use the prolonged systems lZ q + r which are 
obtained by differentiating each equation in lZ q r times totally with respect 
to all independent variables. They are all quasi-linear. If we substitute again 
the power series ansatz into the prolonged system lZ q + r and evaluate at x° , 
we get an inhomogeneous linear system for the coefficients of order q + r. Its 
homogeneous part is determined by the prolonged symbol M. q + r , i.e. the symbol 
of TZ q+r . 

The Taylor coefficients of lower order appear in the matrix and in the 
right hand side of this linear system. Thus we are able to express the coefficients 
of order q + r through the coefficients of lower order. This is the precise meaning 
of constructing a power series order by order. 

This construction will fail, if non-trivial integrability conditions occur, i.e. 
equations of order q + r which are functionally independent of the equations 
contained in the prolonged system lZ q+r and which are satisfied by every solution 
of the system. Such equations arise usually by cross-differentiating and are 
detected only in some higher prolongation. They pose additional conditions on 
the coefficients of order q + r. Hence they must all be known to pursue the 
above described procedure. We call a system which contains all its integrability 
conditions a formally integrable system. 

For formally integrable systems it is thus possible to construct order by or- 
der a formal power series solution. The arbitrariness of the general solution 
is reflected by the dimensions of the prolonged symbols, because at each order 
dim.A/fg+.r coefficients are not determined by the differential equations but can 
be chosen freely [^9|. Formal integrability does, however, not suffice to deter- 
mine these dimensions in advance without explicitly constructing the prolonged 
symbols. This leads to the concept of involution. 

We introduce the class of a multi- index fi = . . . , //„] . It is the smallest k 
for which /i^ is different from zero. If we consider the symbol (0) as a matrix, 
then its columns are labeled by the coordinates . We order them by class, 
i.e. we always take a column with a multi-index of higher class left of one with 
lower class. Then we compute a row echelon form. 

In this solved form the symbol is especially easy to analyze. Since we only 
need linear operations to obtain it, we can always perform the same operations 
with the full system lZ q and thus assume that (||) yields the symbol directly in 
solved form. We denote the number of rows where the leading entry or pivot is 

(k) 

of class k by (3 q and we associate with each such row its multiplicative variables 

X\ , ... j Xfc . 

It is important to note that if we prolong each equation only with respect to 
its multiplicative variables, we obtain independent equations, because each equa- 
tion will have a different leading term. The question is, whether prolongation 
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with respect to the non-multiplicative variables leads to additional independent 
equations. If not we call the symbol involutive. 

Definition 2 The symbol M. q is called involutive, if 

n 

mnkM q+1 = Y. k ^ k) ■ W 

fe=i 

The system 1Z q is called involutive, if it is formally integrable and its symbol is 
involutive. 

(k) 

The above definition of the p q is obviously coordinate dependent. Thus it 
seems, as if the involution of a symbol depends on the chosen coordinate system, 
too. One can, however, show that almost every coordinate system leads to the 

(k) 

same values for the p q . These values are characterized by the property that 
all the sums ^ZILfe f$ , k = 1, . . . ,n, are maximal.^] A coordinate system which 
leads to these values is called 8-regular. Definition || assumes that the f3 q are 
computed in such a coordinate system. Besides there exist alternative methods 
to obtain the correct values intrinsically |^3|, |3l], ^8) . We will return to this point 
in Section [H]. 

The prolongation of an involutive symbol is again involutive. Since prolong- 
ing an equation with respect to one of its multiplicative variables Xi yields an 
equation of class i, we get /3g+ n = J2k=i A? ■ Inductive use of this relation leads 
to 

i=k V / 

and together with Definition to 



rank .M,j+ r = 




Besides the possibility to predict the number of arbitrary Taylor coefficients 
at any order, involutive systems have another advantage compared with formally 
integrable ones. There exists an easily applicable criterion to check whether or 
not a system is involutive. The problem of the definition of formal integrability 
is that one has to prove that a system does not generate non-trivial integrability 
conditions at any prolongation order, i.e. one must check an infinite number of 
conditions. This can, however, be done in a finite manner for systems with an 
involutive symbol. 

Theorem 3 Let lZ q be a q-th order differential equation with an involutive 
symbol A4 q . If no integrability conditions arise during the prolongation of lZ q 
to lZ q+ i, then lZ q is involutive. 

2 Note that this is different from requiring that the /3g fc ' themselves take maximal values! 
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3 Completion to Involution and Arbitrariness 



Since we have seen that involutive systems have many advantages, the question 
naturally arises whether they form only a very special class of systems and 
what to do with a non-involutive system. The interesting answer is given by the 



Cartan-Kuranishi Theorem 18, [23, 28 



Theorem 4 Any system lZ q can be completed to an equivalent involutive one 
by a finite number of prolongations and projections ( i. e. addition of integrability 
conditions) . 

Since this theorem depends on some fairly deep results in the formal theory, 
we will not present a proof but only discuss an algorithm to perform this com- 
pletion. It is based on Theorem || above and consists essentially of two nested 
loops. The inner loop prolongs the system until its symbol becomes involutive. 
The outer loop checks then for integrability conditions and adds them. The 
difficult part of the proof is to show the termination of the inner loop. The 
termination of the outer one follows from a simple Noetherian argument. 

Involution of a symbol can be checked easily using Definition ^|, if we assume 
that the coordinate system is (5-regular what we will do in the sequel. It re- 
quires only linear algebra. Whether or not integrability conditions arise during 
a prolongation, can be deduced from a dimensional argument. 

Denote the projection of the system IZq+i into the q-th. order jet bundle J q £ 
by IZq 1 ^ . Its dimension can be computed indirectly from the identity 

dimT^ 15 = dimftg+i - dimM, +1 (7) 

which reflects the fact that integrability conditions are connected with rank 
defects in the symbol. None has occurred during the prolongation from lZ q 
to lZ q +i, if and only if this dimension is equal to dim 7?^. 

There are essentially two possible reasons for integrability conditions. The 
classical one is that it is possible by some linear combination of equations of 
order q + 1 in TZ q +i to eliminate all derivatives of that order. This is a gener- 
alization of the usual cross-derivative. The other one is that lZ q contains some 
equations of lower order. In order to construct lZ q +i all equations in TZ q must be 
prolonged. If now some equations are of lower order, it might happen that their 
prolongation leads to new independent equations of order less than or equal 
to q. They must be taken into account in the projection to IZ^ . 

Fig. [I] shows this algorithm in a more formal language. 7^+ r denotes here 

(s) 

the system obtained after r + s prolongations and s projections. .M q _{_ r is the 
corresponding symbol. In this form it is comparatively straightforward to im- 
plement it in a computer program. The determination of the dimensions of the 

(s) 

various submanifolds TVq^. r poses the main remaining problem, especially for 
non-linear systems. Ref. [^7], ^8| describe an implementation in the computer 
algebra system AXIOM. 
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[1] 


r <- 0; 




[2] 


compute 72. q +i 


{prolong} 


[3] 


compute M q , Mq+i 


{extract symbols} 


[4] 


until 7Z- q s + r involutive repeat 




[4.1] 


while #multVar(M [ q s l r ) / rankA^' 


l r+1 repeat 


[4.1.1] 


r <— r+1 


{counter for prolongations} 


[4.1.2] 


compute 1Z^l r+l 


{prolong} 


[4.1.3] 


compute M^ r+1 


{extract symbol} 


[4.2] 


if dim7^ r+1 - dim M^l r+1 < dim 


Ktr ^en 


[4.2.1] 


s <- s+1 


{counter for projections} 


[4.2.2] 


compute TZq+ r 


{add integrability conditions} 


[4.2.3] 


compute 1Z q s l r+l 


{prolong} 


[4.2.4] 


compute M q i T ,M i q i r+1 


{extract symbols} 


[5] 


return 7?-S r 





Figure 1: Algorithm for the Cartan-Kuranishi Theorem 



For ordinary differential equations this algorithm becomes very simple. Since 
there is only one independent variable, we find always an involutive symbol and 
cross-derivatives are of course not possible. The only possibility for integrability 
conditions is the prolongation of lower order equations. For partial differential 
equations we recall that the other integrability conditions can always be found 
by considering the prolongations with respect to non-multiplicative variables. 

To conclude this section we briefly recall some results of Ref. j2^] concerning 
the arbitrariness of the general solution which will be needed later. (^) yields 
only the rank of the prolonged symbols, but their dimensions are more inter- 
esting. They can be expressed in a similar way, if we introduce the Cartan 
characters of a differential equation 



^^(^I* -1 )-/^, k = l,...,n. (8) 



They form a descending sequence 



< 5 > 4 2) > • • • > > . (9) 

Now we can write 



dimM g+r = ±afl = ± ( r + k r ^aW. (10) 



fe=i fc=i 



This is the Hilbert polynomial of the differential equation lZ q (it can be written 
explicitly as a polynomial in r). Analyzing the number of arbitrary Taylor 



7 



coefficients in the power series expansion of the general solution and comparing 
with these dimensions yields the following result. 



Theorem 5 The general solution of a first-order system of differential equa- 
tions IZq contains functions depending on k arguments where the numbers fk 
are determined by 

f n = af=m-^\ 
f k = a {k) - a[ k+1) = pf +1) - p[ k) . ( } 

ensures that the fk are always non-negative. Note that Theorem || refers 
to algebraic representations of the general solution, i.e. no integrals or derivatives 
of the arbitrary functions do occur. One can derive more general results covering 
also higher-order equations and more general representations of the solution, but 
we will not need them here. 

We define a gauge symmetry as a fiber-preserving transformation of the 
bundle £ depending on some arbitrary functions of all independent variables 
which maps solutions into solutions. (This implies that f n cannot vanish for 
a system with such a symmetry.) In gauge theories one identifies solutions 
related by a symmetry transformation. In order to obtain information about 
the arbitrariness of the physically relevant part of the solution space we must 
adjust the Cartan characters. 

Let us assume that the gauge transformation can be written in the following 
form 

y ' (12) 
u a = A a (x\u' 3 , \ { °\x), d\2\x), dP\i p) (x)) 

where 70 gauge functions are entering algebraically, 71 gauge functions A« 
are entering through their first derivatives etc. Ref. [|28| shows how one can 
handle more general cases using a pseudogroup approach based on an implicit 
representation of the transformations by differential equations. 

Under this assumption the gauge correction term Aa, ' which must be sub- 

(k) 

tracted from a q to adjust for the effect of the symmetry can be computed 
recursively through 



A < } - + - E few ■ ( 13 ) 

^ ''1=0 i=k+l ^ >' 



where the s k n \q) denote some combinatorial factors, the modified Stirling num- 
bers (earlier called symmetric g-products) introduced in Ref. [p8[ p9|. 
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4 Constrained Dynamics a la Dirac 



Let q % be coordinates in some TV-dimensional configuration space Q. We restrict 
our exhibition to autonomous systems, as explicit time dependencies can always 
be treated by considering the time as additional coordinate in an extended con- 
figuration space. The dynamics of a system is then determined by the condition 
that its action 

S=f L(q\e)dt (14) 

is stationary along trajectories q l (t), where L is the Lagrangian of the system. 
It is well-known from the calculus of variations that this leads to the Euler- 
Lagrange Equations 

d (dL\ dL n 

We pass from the Lagrangian formalism to the Hamiltonian one by a Leg- 
endre transformation. We introduce the canonically conjugate momenta 

* = w (16) 

For regular systems the Legendre transformation provides a one-to-one mapping 
between the velocities q 1 and the momenta pj. In a constrained system this does 
no longer hold; instead one obtains by elimination some primary constraints 

<M<f,ft)=0. (17) 

This implies that not every point of the phase space is accessible for the sys- 
tem (or can be used as initial data) but only a submanifold, i.e. some of the 
coordinates q % do not correspond to true degrees of freedom. 
The canonical Hamiltonian of the system given by 

H c =p i q i -L(q i ,q i ) (18) 

represents no longer the only possible choice. We can add arbitrary combina- 
tions of the constraints without changing its value on trajectories. This leads 
to the total Hamiltonian 

Ht = He + u^(pi (19) 

where the multipliers u e are a priori arbitrary functions of q l ,Pi- 

The constraints must remain stable under the evolution of the system. In- 
troducing the Poisson bracket 

f pM , „. i ,1 dFdG OGdF 
{ W>ft) , G( ,., ft) } = (20) 
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we can express the evolution of any observable F(q % ,pi) concisely 

F = {F 1 H T }. (21) 

Thus we are lead to the requirement 

W/,5r}«0. (22) 

The sign rs signals that this is a so-called weak equality, it must hold only af- 
ter taking all constraints into account. By a standard argument in differential 
geometry this implies that the Poisson bracket in (^2|) must be a linear com- 
bination of the constraints. There are three possibilities for (p2|): (i) it yields 
modulo the constraints an equation of the form 1 = 0; (ii) it becomes = 0; 
(hi) we obtain a new equation ip(q l ,pi) — 0. 

(i) means that our equations of motion are inconsistent. This implies that 
they do not possess any solution. Hence the Lagrangian is physically invalid, 
(ii) is of course the desired outcome, (iii) results in a secondary constraint. It is 
added to the other ones. We must of course then check whether all secondary 
constraints remain stable under the evolution of the system, i.e. we have to 
repeat the procedure until we either encounter case (i) or all constraints lead to 
case (ii). This is the so-called Dirac algorithm. 

If secondary or higher constraints occur, we must distinguish whether or 
not they depend on the multipliers xfi. If yes, we can solve for some of them 
which are then no longer arbitrary. This indicates the presence of second class 
constraints, as a first class constraint tp Poisson commutes weakly with all other 
constraints fa, i.e. 

{V>,M~0. (23) 

It is well-known that first class constraints generate gauge symmetries. Sec- 
ond class constraints correspond to unphysical degrees of freedom; a typical 
example is the pair q 1 — and p\ = 0. These unwanted degrees of freedom can 
be eliminated using the Dirac bracket. Let Xi denote all second class constraints 
and define the matrix C by 

(\. (24) 
This matrix is always non-singular and we can define 

{f,9\ = {f,9}-{f,Xi}{C- l r{X3,9} ■ (25) 

In the canonical quantization of the system the Dirac brackets and not the 
Poisson brackets are transformed into commutation relations. 

One of the fundamental goals in constrained dynamics is to count the number 
of degrees of freedom of the system. If there are Np first and N$ second class 
constraints in the system, then the number F of dynamical degrees of freedom 
is given by 

F = N-N F - ^N s . (26) 
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This simply reflects the fact that two second class constraints are necessary to 
eliminate a degree of freedom, as we need one for the coordinate and one for 
the momentum. A first class constraint leads to a symmetry and thus to an 
arbitrariness in a coordinate (or a momentum). A real elimination requires a 
gauge fixing condition, i.e. we add a new constraint which turns the first class 
constraint into a second class one. 



5 Involution Analysis 

Now we will analyze the equation of motions from the point of view of formal 
theory. Let the bundle £ be given by £ — Q x T, where Q is the TV-dimensional 
configuration space with coordinates q % and T the time axis, together with 
the natural projection £ — * T. ( [i"5| ) represents a second-order equation whose 
symbol is determined by the Hessian matrix 



d 2 L 

M « = WW (27) 



If the symbol has rank N, the Euler-Lagrange Equations are normal and no 
constraints occur. Its general solution is parameterized by 2N arbitrary con- 
stants. 

If, however, the symbol has lower rank, it is possible to eliminate the second- 
order derivatives in some of the equations. Now it is no longer obvious, whether 
or not ( |i"5| ) is involutive. Since we are dealing with an ordinary differential 
equation, the symbol is always involutive. But the prolongation of the obtained 
differential equations of lower order might lead to integrability conditions, if the 
resulting equations are independent of the remaining second-order equations 
in ([15]). Then we have to check whether some of these conditions are again of 
lower order; in that case we have to repeat the procedure. 

After a finite number of iterations we will obtain either an inconsistency or 

(s) 

an involutive system TZ^ of the following form 

q j =f j (q i ,q i ,q n ), j = 1, . . . , J, n = J + 1, . . . , N , 
K<2 ] : { q k =g k {q\q n ), k = l,...,K, n = K + 1, . . . , N , (28) 

j m = h m (q n ), m = l,...,M, n = M + 1, . . . , N , 

with M < K < J. Refs. |33[ |34) contain detailed treatments of constrained 
systems in the Lagrangian formalism. A closer look reveals at once that it 
corresponds exactly to the completion algorithm presented in Section |3] applied 
to a system of second-order ordinary differential equations. Zeroth and first 
order equations are called there constraints of A and B type, respectively. 

To relate our approach to the standard one by Dirac we pass again by a 
Legendre transformation to the Hamiltonian formulation. At the level of the 
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differential equations this means that instead of the configuration space Q the 
phase space P is used to construct £. In other words, we introduce N additional 
dependent variables p, and transform ( |l5| ) into the first-order equation 

- f*L 

7? j : { d $ L i = l,...,N. (29) 

= w 

The first set of equations consists of course just of the Euler-Lagrange Equations 
with the second set of equations plugged in. 

If the matrix Mjj has full rank, the second set of equations can be solved 
for the q 1 and (^9|) is a normal equation. Otherwise we obtain some algebraic 
equations of the form 4>i{q l ,Pi) = 0. They are of course the primary constraints 
of Dirac. Just as in the Lagrangian formulation, involution of ( p9| ) depends on 
the behavior of the prolonged equations 

Dt ^ = W q + Wi Pi= (so) 

We are only interested in these equations restricted to IZi . This corresponds to 
the weak equalities used in the last section. Like in the Dirac algorithm there 
are three possibilities for the result of the restriction of each of the equations 
in (|30|): (i) it yields an inconsistency; (ii) it vanishes identically; (iii) we obtain 
a new independent equation. 

If secondary constraints appear, we must repeat the procedure to check the 
consistency of the equations of motion (|29|). After a finite number of steps we 
will cither have found an inconsistency or we will have constructed an involutivc 
equation of the form 

dL ■ 1 

Pi = g-j, l=l,...,N , 

K[ S) : \ gj = fj( q *,q", Pi ), j = 1, . . . , M , n = M + 1, . . . , N , ( 31 ) 
fa{q i ,p i )=Q, £=1,...,K. 

In the Dirac approach one does not use the Hamiltonian equations (p9|), 
but one introduces some multipliers and takes those derived from the total 
Hamiltonian To justify this we look at the differential of the canonical 

Hamiltonian (Hq) 

dH c 

Thus on Hi we obtain 



(32) 



dL 

dHc\ ni ^q l d Pl ~—dq\ (33) 
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Two one-forms which coincide when restricted to the constraint surface, i.e. 
the submanifold of £ defined by the constraints 4>e = 0, can differ only by a 
linear combination of the form u e d(j)z with arbitrary coefficients u . Since 



dHr 



dH c 



dq 1 



dH c 

dpi 



dpi , 



(34) 



we obtain the following equations of motion (sometimes called Hamilton-Dirac 
equations) living in an extended phase space 



dH, 



c 



Tlx ■ { 



-Pi = 



dpi 
8H C 



dpi ' 



dq 1 dq 1 



(35) 



Here the coefficients u l must be considered as additional functions of t or in 
the language of differential equations as additional dependent variables. ( ^9|) is 
obtained, if we use the first set of equations in ( fjB]) to express u l through q l , q l . 
Thus both systems are equivalent. 

The Dirac algorithm is equivalent to the completion of system (|35|). It 
requires the analysis of the prolongations of the constraints restricted to TZ\. 
They can be concisely written using Poisson brackets 



D t (, 



OHc 



dH 



c 



dq 1 



dpi 
H c } 



dpt dq 1 
u k $y4>u4> k } 



b e d<j) k d<j>i d<f) k 



dq 1 dpi dpi 



(36) 



Actually, to obtain a full equivalence we should write the multiplier as deriva- 
tives u l — v l . This is only important in the case that second class constraints 
are present. Then we obtain equations determining some of the u e . In principle 
we must prolong then these equations, too. This unnecessary step which yields 
no new information can be omitted by using derivatives. This appears also from 
a physical point of view somewhat more natural, as the multipliers correspond 
to velocities. 



6 Counting Degrees of Freedom 

The classical expression ( p6| ) for the number of degrees of freedom depends on 
the distinction into first and second class constraints. This requires, however, 
the introduction of a Poisson structure. We will show now that it is possible 
to obtain an intrinsic expression for this number without performing such a 
distinction. 
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We will start with the Hamiltonian equations of motion (|T]). In order to 
obtain the number of degrees of freedom we must count the number of constants 
necessary to characterize a physical state. Since dim J\£ = AN, we obtain 

(s) 

dim 7^ = 3N —M —K . Thus in a power series expansion of the general solution 
this number of zeroth and first order coefficients can be chosen arbitrarily. 

If we identify the arbitrary functions with q l for i > M, as these are not 
restricted by (|3|), we must subtract N — M constants coming from the arbitrary 
functions. Thus an initial state at t = to is specified by 2N — K constants. 
Depending on the choice of the arbitrary functions we will, however, obtain 
different values for q l {t\) andpj(ti) at some later instant t\. As these correspond 
nevertheless, by definition, to the same physical state, we must subtract further 
N — M constants for the gauge symmetry. Thus a physical state is specified by 
N + M — K constants and the number of degrees of freedom F is half of this 
number 

F=-(N + M-K). (37) 

We find for u[ s) that /3[ 1] = N + M and hence a[ 1] = N -M. Expressing 
M by N and and similarly K by N, and dim 72.^ yields an intrinsic 

( s) 

expression for F independent of any specific representation of the manifold 1Z\ 

F =idim^ s) -af ■ (38) 

If we use instead of the Hamiltonian equations of motion (|2^) the Euler- 
Lagrange Equations (|l5|), we obtain the analogue result 

F = - dim74 s) - 4 X) . (39) 

Both expressions yield always the same value, as we will obtain exactly the 
same dimension and Cartan character for the final involutive equation no matter 
whether we work in first- or second-order because of the different dimensions of 



the base spaces |23|, |28 . 

As a by-product this implies that J in (|||) equals M in ( |3l| ) and the sum of 



K and M in ( |28[) equals K in (|3l|). Thus we have the same number of constraints 
in both approaches, if we omit the introduction of multipliers. This might not 
be very surprising, we will, however, see later that this does no longer hold in 
field theories. 

We hope to study in a future paper the distinction into first and second class 
constraints in more detail. Then we will also be able to discuss the relation 
between these results and the classical formula (|26|). For the moment we just 
note that according to Theorem || the general solution contains arbitrary 
functions. In the classical terminology this arbitrariness stems from the gauge 
transformations generated by the primary first class constraints |l4], [ll| . Thus 
their number is a^. 
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7 Example 



Consider the classical problem of a particle whose movement is restricted to the 
surface of a sphere in a D-dimensional space but otherwise free, often also called 
the rigid rotator ||. Without loss of generality we can take the radius of the 
sphere as one and start with the Lagrangian 

L(q\q\\,\) = jq 2 + \(q 2 -l). (40) 

(q 2 = q l qi, etc.) A is here obviously a multiplier. The canonically conjugate 
momenta pt , it are given by 

Pi=mqi, i = l,...,D, ^ 

7T = . 

If we introduce a further multiplier /z, we can write the total Hamiltonian as 

Ht = - Ml 2 - 1) + ^ (42) 

Obviously there is one primary constraint, namely 7r = 0. The next three 
steps of the Dirac algorithm lead to the constraints q 2 = 1, pq = 0, and finally 
p 2 = —2mA. It is easy to see that all these constraints are second class. The 
system contains thus D — 1 degrees of freedom which can be formally calculated 
by subtracting from the dimension of the configuration space — D coordinates qi 
plus one coordinate A — half the number of second class constraints, i.e. 2. 

Now we will obtain the number of degrees of freedom using a formal analysis 
of the Euler-Lagrange Equations 

r, : ; ;0; (43) 



The completion to an involutive equation requires four projections leading to 
the integrability conditions qq = 0, qq + q 2 = 0, A = and finally A = 0. After 
some trivial manipulations we have thus the involutive equation 

mq\ - 2\q t = , A = , 
^ 4) : \ qq = 0, A = 0, (44) 

q 2 - 1 = , mq 2 + 2A = . 




It is easy to see that this represents a finite type equation and thus there are no 
first class constraints. Since dim = 2D — 2, ( |3^ ) yields F = D — 1 degrees 
of freedom, in perfect agreement with the Dirac analysis. 
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Alternatively, we can analyze the Hamiltonian equations of motion 

1=0, 



Hi 



Pi - 2Xqi = , 
mqi - pi = 0, 



7r — q- 

7T = 0. 



(45) 



Again the system becomes involutive after four projections with integrability 
conditions 7r = 0, qp — 0, p 2 + 2mA = 0, and finally A = 0. This yields 



n 



(4) 



pi - 2Xqi = , 
mqi -pi = 0, 
TT = 0, 
qp = 0, 



7T = 0, 

A = 0, 



P 



2mA = , 
1 = 0. 



(46) 



The analysis of this equation leads of course to exactly the same number of 



degrees of freedom, as dim TZ\ — dim IZ^' and both are finite type equations 



(4) 



8 Field Theories 

In the section on point mechanics we studied three different ways to write the 
equations of motion: the Lagrangian equations, the Hamiltonian equations ob- 
tained from the latter one by direct application of the Legendre transformation, 
and finally the Hamilton-Dirac equations where one includes multipliers for the 
primary constraints. In field theories there is even more choice, as there exist 
at least two different ways to perform the Legendre transformation. 

The standard way entails the explicit choice of a time variables and leads to 
a non-covariant formalism. From the point of view of differential equations this 
approach has a further disadvantage: it is not a truly first-order formalism, as 
the Hamiltonian field equations will generally still contain second-order spatial 
derivatives. Thus for the application of involution theory it is probably more 
appropriate to use the so-called De Donder-Weyl approach Jl6[ |25) which leads 
to a covariant first-order formalism. 

We will therefore restrict ourselves in the sequel to the analysis of field 
equations in Lagrangian formulation and leave the discussion of the Hamiltonian 
approach for the future. This suffices for the purpose of this article. 

Many articles on the theory of systems with constraints have the following 
structure: The theoretical results are derived in the finite-dimensional case, i.e. 
in point mechanics; the examples and applications stem, however, from field 
theories. The connection is made with a remark like "The generalization of 
these results to field theories is straightforward" . But this point of view is a 
bit optimistic, as a more careful discussion (see e.g. Q) reveals many subtle 
problems. 

Although on the surface the main difference lies in the fact that Poisson 
brackets are now computed via functional derivatives instead of partial ones, 
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many elementary concepts in the finite-dimensional theory become rather tricky 
in an infinite-dimensional setting. For instance linear combinations must now 
be substituted by integrations. But to require that an integral vanishes is a 
much less stringent condition than the vanishing of an algebraic expression and 
depends decisively on the considered function space. 

Similarly, inverses as they are used in the construction of Dirac brackets are 
no longer uniquely defined. Often already the distinction into first and second 
class constraints can be rather problematic and statements like the number of 
second class constraints is always even do not make sense any more. One must 
introduce the new concept of proper and improper constraints |IJ . Further prob- 
lems stemming from the choice of coordinates will be discussed in Section [k| 

To really solve these problem one must usually resort to fairly complicated 
methods from functional analysis. We will concentrate in this section, however, 
on another point: The naive generalization of the Dirac analysis does not cor- 
respond to the completion to involution of the field equations. Thus in general 
it does not suffice to prove their consistency. 

Since field equations are partial differential equations, involution becomes a 
more complicated concept. The prolongation of lower order equations represents 
no longer the only way to generate integrability conditions. Constraints are 
mostly equations of lower class. In a typical field theory the base space X 
of independent variables is a D-dimensional space-time. We can identify the 
variable x D with time. Thus in the usual terminology, the equations of class D 
are the evolutionary ones; the remaining ones are constraints. 

The naive generalization of the Dirac algorithm prolongs all constraints only 
with respect to time, as it relies solely on Poisson brackets with the Hamiltonian. 
If all constraints are of class D — 1, this corresponds to our approach, because 
it suffices to analyze the prolongations with respect to the non-multiplicative 
variables and we find in this case only one, namely x . 

But now the question arises as to what happens if constraints of lower class 
appear. Then prolongations with respect to the other multiplicative variables, 
i.e. with respect to some spatial coordinates, may lead to additional integrability 
conditions not considered by this naive approach. 

In order to exhibit this effect in "pure form" we begin with a highly unphys- 
ical example without any kinetic term in the Lagrangian density. But we will 
later show that this is not the important point. Consider the class of systems 
described by the Lagrangian density 



on a three-dimensional flat space-time with coordinates x,y,t. The fields A, /i 
represent again multipliers, whereas /, g denote fixed but so far arbitrary func- 
tions. Variation with respect to </> yields the equation 



A, fx] = »[d x cj> - /(</»)] + X[d v <P - g(</>)] 



(47) 



d xf i + d y \ + f((l>)v + g'{<t>)A = 0. 



(48) 
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More interesting equations are obtained from the multipliers. They generate an 
over-determined system for <j) 



dx4>-f{4>) = 0, 

V y ' (49) 

d v 4>-g{<t>) = Q. 

Obviously this system is consistent, if and only if f.g satisfy the equation 

f'g = fg'- (50) 

This requires that / is a multiple of g. It is easy to see that only under this 
condition the Euler-Lagrange Equations are involutive. We must conclude that 
most of the Lagrangian densities (|4^) are physically invalid. 

Now we look at the outcome of the naive Dirac algorithm applied to this field 
theory. Obviously all three canonically conjugate momenta 71"^, 7r^, tt\ vanish and 
represent thus primary constraints. The total Hamiltonian density is given by 

H T = -fJt[d x (p - f((j>)} ~ \[dy<t> - g{(f>)) + + VTT^ + W7T\ . (51) 

This yields the following secondary constraints 

{iC4,H T } = ~[d xf i + d y X + f'Wn + g'(<p)\] , 

K,H T } = 9^-/(0), (52) 

{tta,Wt} = d v (j)-g{4>). 

Thus we obtain exactly the Euler-Lagrange Equations above. But note the 
crucial difference in the further analysis. Following the naive Dirac analysis 
we look only whether these equations are algebraically related, i.e. whether one 
vanishes, if we take the others into account. But this does not happen here. 

Since all constraints are second class, we continue to compute the tertiary 
constraints in order to fix the multipliers u, v, w introduced in the total Hamil- 
tonian density. They are 

3 x v + d y w + f"(4>)fm + f{<t>)v + g"{ct>)\u + g'{4>)w = , 

d x u- f'(cp)u = 0, (53) 

d y u — g'{4>) u = . 

We find the interesting phenomenon that although there are six second class 
constraints for three degrees of freedom A, they do not fix all multipliers. 
v, w appear only in the first equation, hence one of them can be chosen arbi- 
trarily, u is the solution of an over-determined system which is consistent, if / 
and g satisfy (^0|) . 
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This example shows that the real failure of this approach lies in the purely 
algebraic treatment of the constraints. In point mechanics we could always 
substitute one constraint in another one, if both contained the same coordinate 
q l or a (time) derivative of it. This is no longer possible in field theories, 
because there might be derivatives with respect to other coordinates present. 
To check the consistency of such constraints requires the analysis of integrability 
conditions. 

This analysis is surely trivial in a primitive example as above. But in more 
complicated cases it is rather difficult to decide when one has found all inte- 
grability conditions. Assume for instance that a field theory in 3+1 dimensions 
leads (among others) to the following constraints for some field <j) 

4>zz + wt>xx = , 

<j)yy=0. 

In this famous example due to Janet one needs five prolongations and two 
projections to show that there are exactly two integrability conditions, namely 
4>xxy = 4>xxxx = 0. Such problems lead at the end of the last century to the first 
steps towards the development of the formal theory! 

One might argue that ( fl7|) is a rather peculiar Lagrangian density. But the 
addition of a kinetic term for <f> does not really change the outcome, although 
the computations are slightly more complicated due to the appearance of further 
integrability conditions. Take 

C = l(d t & 2 +£. (55) 
Its Euler-Lagrange Equations are 

d tt <t> + d xf x + d y X + /'(»// + g'(<t>)\ = , 

dx.4> - f{4>) = o , 

d y <P-g(<P) = 0. 

Only the first equation has changed. But since the only second-order deriva- 
tive involves <j>, our completion procedure generates besides ( |50| ) two further in- 
tegrability conditions involving second-order derivatives of /x, A. We omit them 
here, because they are rather complicated. 

We also do not show the Dirac analysis of this system. It suffices to note 
that the two new conditions (in Hamiltonian form) are also found by the Dirac 
algorithm. But this was to be expected, because their construction requires 
prolongations with respect to time. Only ( jspf ) is again overlooked, for no spatial 
cross-derivatives are performed. 

The natural question is whether systems of this form are unphysical for some 
reason or whether this effect occurs often and might lead to a wrong results. 
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(54) 



(56) 



Such purely spatial cross-derivatives (or to be more precise: linear combinations 
of spatial prolongations) are only necessary, if the field equations form an over- 
determined system. This might appear at first sight quite unusual but one can 
find such systems in the literature (see e.g. Q). A typical example is a gauge 
fixing condition built into the Lagrangian density. 



9 Higher- Order Lagrangians 

Theories described by Lagrangian containing higher order derivatives |], [D| 
p0| , pi[ can be treated in exactly the same way. Using again the notation of 
Section ||such a Lagrangian is a function L(x z ,u a ,p'^). The x l are coordinates 
on the underlying space-time in the case of a field theory or just the time in 
point mechanics. Similarly, the u a denote either the fields in the theory or the 
generalized coordinates and the represent the derivatives. 

From the calculus of variations it is well-known that the Euler-Lagrange 
Equations can now be written using the Euler operators ]2jJ 



If /i = [jtii, ... , /ij, then = D^ 1 ■ ■ ■ D£ n where Di denotes the total derivative 
with respect to x l . Obviously we recover ( |l5"|) , if there is only one x % , namely 
the time t, and L depends only on derivatives with = 1. The sum in (^) is 
always finite, as L contains only derivatives up to a given order. 

As in the standard case there exist at least three possibilities for the starting 
point of the involution analysis. The simplest choice is of course to use directly 
the Euler-Lagrange Equations (|57|) . Alternatively one can pass to a Hamiltonian 
formulation. This requires now the introduction of several momenta conjugate 
to each field. Then one can either transform the Euler-Lagrange Equations 
directly or one can derive in the same manner as before the Hamilton-Dirac 
Equations by introducing multipliers. 

There is no need to repeat the arguments of Section ^|, as they still apply 
in the same way. Instead we consider as an example Podolsky's generalized 
electrodynamics in the Lagrangian formalism. It demonstrates the typical 
way an involution proof proceeds for field theories. Refs. [Q |3(| contain further 
examples like Yang-Mills or Einstein Equations. The Lagrangian density is given 

by 

C = ~\f^ v - a 2 d p F°Pd T F; . (58) 

The space-time indices run from 1 to D; we identify x D with the time, a is a 
constant. If it vanishes, we recover the standard Maxwell theory. 
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In terms of the vector potential A p the Euler-Lagrange Equations ( p7[ ) give 
a fourth-order system 

Ki : (1 - 2a 2 u)uA ll - p [(1 - 2a 2 u)d v A v ] = 0, n = l,...,D (59) 

with the D'Alembertian □ = rj^ v d p d v . 

According to Definition || we must first check whether the symbol M.^ is 
involutive. It is given by the equations 

M, : rrrrv%,p* ~ V° T V a(3 v 5 OT0S = , a = 1, . . . , D (60) 

where V^ upa is a placeholder for the derivative d^ pa A a . It is easy to see that 
for a 7^ D we can choose in each equation the variable vf) DDD as pivot, i.e. all 
these equations are of class D. For a — D, however, the corresponding variable 
cancels. We can obtain at most a pivot of class D — 1, e.g. vjj_ 1 D _ 1 D _ 1 D _ 1 . 
Thus 

Pi D) =D-l, /3f- X) =l, /3f=Q. (61) 

One can prove that this cancellation for a = D does not simply stem from 
a singular coordinate system either by using the tableau of the system (cf . [£8| ) 
or by arguing that = D is not possible, as the system has a gauge symme- 
try |n| . This argument relies on the results of Ref. |2£| on the arbitrariness of 
the general solution. 

In order to apply Definition ^| we must compute next the rank of the pro- 
longed symbol A4$. It is defined by 

M 5 : ^V CT ^ P , 7 -»7 OT ^ o % 75 = 0, a, 7 = l,...,D. (62) 

It is easy to see that these equations are not all independent, because if we set 
a = 7 and sum the result vanishes. 

It follows from the discussion before Definition and the obtained values for 

(k) 9 

the PI ' that rank_M 5 > D z - 1. But since we have already found an identity, 
the same value represents an upper bound for the rank. Thus the rank must be 
exactly D 2 — 1 and by (f|) Ma is involutive. 

In order to see whether integrability conditions occur we must check whether 
this identity holds only at the level of the prolonged symbol or also if we use 
the full prolonged equations. We know already that the fifth order derivatives 
cancel, but it could happen that some lower order equation remains. It is, 
however, easy to see that this is not the case. This is analogue to the Noether 
identity in the Maxwell theory. Hence IZa is involutive. This implies that in 
the Lagrangian formalism only one constraint appears, namely the equation for 
a = D which is of lower class. 

For later use we note the Cartan characters of the theory 

4 4) = i , 4 3) = 15 1 4 2) = 40 , 4 1} = 80 • (63) 
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Adjusting for the symmetry — > A M +9 M A with (|T^) yields the following gauge 
corrected values (71 = 1) 

af = , 4 3) = 10 j 4 2) = 25 1 "4 1} = 45 • ( 64 ) 

The Hamiltonian treatment is similar but more involved, as now secondary 
constraints appear. There is no need to detail it here. Since the only constraint 
in the Lagrangian formalism is of class D — 1, the Dirac analysis is sufficient 
and equivalent to the involution analysis. Thus we just recover the calculations 
presented in Ref. |jio| . 

The constraint analysis of Podolsky's generalized theory is very similar to 
the standard Maxwell theory. In both cases we find that the Euler-Lagrange 
Equations are already involutive due to the Noether identity, whereas in the 
Hamiltonian formalism we must perform a few steps until we reach an involutive 
system. The same effect can be observed in other field theories. This seems to 
imply that at least at classical level the Lagrangian formalism is more efficient, 
as it yields faster an involutive system. 



10 "Field Theoretical" Degrees of Freedom 

The classical procedure to count degrees of freedom in field theories is simply 



to stick to the rule (26) used in the finite-dimensional case. N denotes now 
the number of fields. The argument is essentially the same: Each constraint 
"fixes" one field in the phase space and in the case of a first class constraint the 
symmetry eliminates a further degree of freedom. 

The problem with this approach is that in field theories constraints are usu- 
ally differential equations. Hence they cannot really fix a field; there remains 
some freedom. Speaking about arbitrary functions in the general solution the 
idea behind the rule above seems to be that a degree of freedom corresponds to 
the possibility to prescribe as initial data a function of D — 1 arguments, i.e. 
of all spatial variables. This is for instance the case in a regular theory, as its 
field equations form a normal system satisfying the conditions of the Cauchy- 
Kowalevsky Theorem J9] . 

A generalization of this idea was already introduced by Einstein with his 
definition of the strength of a system of differential equations. Up to a numerical 
factor depending on the dimension of space-time it can be identified with the 
number of arbitrary functions of D — 1 arguments [ p9| provided the system 
is absolutely compatible. Latter condition entails that there are no arbitrary 
functions of D arguments. 

Of course such considerations make sense only after taking gauge symme- 
tries into account, because such a symmetry leads to arbitrary functions of all 
independent variables. Einstein's definition of the strength contains such cor- 
rection terms. Based on the formal theory of partial differential equations we 
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propose as an intrinsic definition for the number of "field theoretical" degrees 
of freedom the gauge corrected Cartan character or q D . This definition covers 
also theories derived from Lagrangians containing higher-order derivatives. 

According to Theorem || a q D corresponds to the number of arbitrary 
functions of D— 1 arguments in the general solution modulo the gauge symmetry. 
Of course we assume here that a q — as in Einstein's approach. Otherwise 
either the field equations are under-determined or the gauge group was not 
correctly identified. For a regular theory no gauge correction is needed. In this 
case a q — and a q D — m where m is the number of fields. Thus we recover 
the usual result. 

For the Maxwell theory one obtains the following gauge corrected Cartan 
characters Q 

4 4 >=0, 4 3) = 4 , 4 2) =6, 4 1)=0 - (65) 

(The same values can be obtained without gauge correction by directly analyzing 
the field equations in field strength formulation.) Thus we obtain in perfect 
agreement with the usual result 4 degrees of freedom. In contrast Podolsky's 
generalized theory possesses 10 degrees of freedom as evident from (|64|). 

It is important here to note that the Cartan characters are intrinsically 
defined and thus independent of any specific coordinate system. This is of 
course a property one should expect of a reasonable definition for the number 
of degrees of freedom. As Stcinhardt Q showed with some explicit examples 
(see also the discussion in |[4]), the classical approach encounters problems, if 
"wrong" coordinates are used. 

He considered among others the simple example of a free massive scalar 
field in 1+1 dimensions. In standard coordinates this system described by the 
Lagrangian density 2C = —d^tpd^cf) — m 2 4> 2 is obviously regular and contains 
one degree of freedom. In light-cone coordinates x^ = (x±t)/2 the Lagrangian 
density becomes 

C = d + (j)d-cl>-m 2 4> 2 l2. (66) 

If we choose x + as new evolution parameter, the canonically conjugate momen- 
tum is 7r = d^4> and independent of the velocity d + 4>. Hence the system is 
constrained. 

It is quite subtle to decide whether this constraint is first or second class, 
but here we are not concerned with these difficulties. The important fact is that 
there appears a constraint and hence according to ( p6| ) the system has less than 
one degree of freedom! Obviously this cannot be correct. We are not aware of 
any proposal in the literature for a modified formula to count degrees of freedom 
that takes this effect into account. 

In the involution analysis a similar phenomenon occurs. As already Stein- 
hardt pointed out, the appearance of the constraint is intimately connected with 
the fact that the light-cone coordinates are the characteristics of the field equa- 
tions. Actually this represents a special case of a more general problem, namely 
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the (^-regularity of a coordinate system fl23| which was already mentioned in 
Section 0. 

If a coordinate system is not <5-regular, the procedure to compute the (3 q 
described in Section || yields too small values. This corresponds to the too many 
constraints found by Steinhardt. For instance in our case the field equations are 

2d+d-(j) + m 2 (t) = 0. (67) 

Obviously there is no derivative of class 2 in this equation. But it can be 
generated with a simple coordinate transformation, namely going back to the 

(2) 

original ones. Thus the correct value for is one and not zero. This yields 
= and = 1- Since no gauge correction is necessary here, we obtain 
the expected result: one degree of freedom. 

There exists a simple method to determine the correct values of the (3^ in 
any coordinate system without performing a coordinate change. It makes use 
of the generalized tableaux of a differential equation. Their ranks provide an 

(k) 

intrinsic definition for the p q ; i.e. their determination is a simple problem in 
linear algebra. For lack of space we cannot detail this approach here but refer 
to the literature |5l], ^8|. The important point is that our definition for the 
number of degrees of freedom can be applied to any system in any coordinate 
system and leads always to the same number. 



11 Conclusion 

Obviously it is one of the most elementary requirements on a system of differ- 
ential equations to be consistent, i.e. to possess at least a formal power series 
solution. The formal theory of differential equations provides us with a powerful 
tool to check this property: the involution analysis. We have shown how it can 
be used in a physical context, namely in constrained dynamics. Here the dif- 
ferential equations arise typically as the equation of motions derived from some 
Lagrangian L. 

One should keep in mind that the motivation for the Dirac analysis is ex- 
actly the same. Its first task is less to exhibit all constraints but to prove the 
consistency of the equations of motion. Hence it is not very surprising that we 
find that for finite-dimensional systems the Dirac analysis coincides with the in- 
volution analysis. Both approaches yield (in some sense merely as a by-product) 
all constraints or integrability conditions, respectively, of the system. 

Although it seems to be a commonly accepted claim that the Dirac analysis 
can be extended without modifications to field theories, we have given examples 
where in our opinion this classical approach is not sufficient. Their construction 
was based on the simple observation that Dirac takes only the temporal evo- 
lution of the system into account. He does not consider spatial prolongations. 
Hence his approach is incomplete and not able to prove the consistency of the 
field equations. 
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One cannot really speak of a failure of the Dirac approach, but one must note 
that it must be augmented by some kind of analysis of the spatial integrability 
conditions. One could for instance use an approach like the involution analysis 
for this purpose. But we think that it is conceptually easier to use one method 
instead of a combination of several ones. 

We believe that the involution approach is more flexible than the Dirac 
method. It can be applied to the Lagrangian as well as to the Hamiltonian 
formalism. In a future article we will show that the so-called symplectic method 
for first-order Lagrangian js| also fits into this scheme. The same holds for 
higher-order Lagrangians. As soon as the equations of motion are obtained, 
either by a variational principle or some other way, the involution analysis can 
start. This is especially important for systems with anholonomic constraints 
which cannot be treated with the usual variational methods pf| . 

An important advantage of the involution analysis is that it represents a 
geometric framework, i.e. all constructions are intrinsic and coordinate inde- 
pendent. This is not true for the Dirac approach. It is well-known that the 
number of constraints can be different in different coordinate systems. A typi- 
cal example are light-cone coordinates. In contrast, our definition of the number 
of degrees of freedom for a field theory via the Cartan characters is completely 
intrinsic. To our knowledge no such definition has so far been proposed in the 
literature. 

One might wonder where all the subtleties of the Dirac analysis like reg- 
ularity conditions on the constraints, ineffective or reducible constraints have 
disappeared p5[ . They are of course still present. Most of them are hidden 

— (s) 

behind the calculation of the dimension of the various submanifolds TZ q + r used 
in the completion algorithm depicted in Fig. [I]. But this is a classical problem 
in geometry and for special types of constraints there may exist alternative ap- 
proaches. For instance polynomial non-linearities represent probably the most 
important case in applications. For them we can avoid a discussion of most of 
the mentioned problems by using Grobner bases techniques |]l7|] . 

Another effect which might lead to problems is that the rank of the symbol 
(or more generally the numbers l3q k ^) does not need to be constant. It might 
change, if certain additional differential equations hold. This generalizes the 
classification of the eigenvalues of the Hessian of the Lagrangian as given by 
Lusanna [fl9| . 

Finally, one should note that by making contact with the formal theory of 
differential equations one obtains suddenly a well understood object, namely 
an involutive system. Many properties of such systems are known and many 
techniques have been developed for their further analysis. All these results are 
now available for constrained systems. 
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